Conformal Invariance and SLE in Two-Dimensional Ising Spin Glasses 
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We present numerical evidence that the techniques of conformal field theory might be applicable 
to two-dimensional Ising spin glasses with Gaussian bond distributions. It is shown that certain 
domain wall distributions in one geometry can be related to that in a second geometry by a conformal 
transformation. We also present direct evidence that the domain walls are stochastic Loewner (SLE) 
processes with « w 2.1. An argument is given that their fractal dimension df is related to their 
interface energy exponent 9 by df — 1 = 3/[4(3 + #)], which is consistent with the commonly quoted 
values df w 1.27 and 6 w -0.28. 

PACS numbers: 11.25.Hf, 75.50.Lk, 05.50.+q 



The powerful tools provided by conformal field theory 
(CFT) have permitted the determination of the expo- 
nents associated with most two-dimensional (2d) critical 
phenomena Q . Unfortunately, CFT has not to date pro- 
vided any results on systems like spin glasses, which in 
two dimensions have a zero temperature transition i.e. 
T c = 0. As the temperature of the system T is reduced 
to zero, the correlation length £ increases to infinity 0,0 
as £(T) ~ 1/T V '. The Hamiltonian of the system is 
H = —J2<ij> JijSiSj, with S% = ±1. If the nearest- 
neighbour bond distribution J^j is continuous, as in the 
Gaussian distribution, then the exponent v = —1/9, 
where 9 is the exponent which describes how the energy 
AE of a domain wall (DW) which crosses a system of lin- 
ear extent L depends upon L: AE ~ L 6 3]. According 
to numerical studies 0,0 9 lies between -0.28 and -0.29. 
The DW has a fractal dimension d f « 1.27 ± 0.01 @. 
While the disorder present in spin glasses means that at 
a local level there is not even translational invariance, the 
existence of a diverging lengthscale such as £ (T) suggests 
that on long lengthscales such microscopic features could 
become irrelevant and that possibly even CI might arise. 
It is the chief purpose here to provide numerical evidence 
that this is indeed the case and so hopefully pave the way 
to eventually determining exponents like 6 using CFT. 

In the first part of this letter we present a numerical 
study of whether there is conformal invariance of the DW 
distribution in 2d spin glasses. Within our numerical ac- 
curacy, CI does seem to hold in the thermodynamic limit. 
Then we next present numerical evidence that DWs in 2d 
are stochastic Loewner evolution (SLE) processes 0, and 
finally we suggest a relationship between 9 and df. 

CI of the DW distribution implies that given two ge- 
ometries related by a conformal transformation, then the 
probability of finding the DW in a given configuration in 
one geometry is related to the probability of finding the 
conformally transformed DW configuration in the trans- 



formed geometry. We will find a transformation F(z) 
mapping geometry (a) onto (d) in Fig. ^ The two rect- 
angles have periodic boundary conditions so that the left 
and right edges are identified, while the top and bottom 
edges are open i.e. they have the topology of an annulus. 
The two slits in (of) also have open boundary conditions, 
so that no bonds cross the slits. The rectangle (a) has 
an arbitrary aspect ratio; we can tune the distance be- 
tween the slits in (d) as desired, with the aspect ratio of 
the right rectangle being a function of the distance be- 
tween the slits and the aspect ratio of the left rectangle 
as given implicitly by Eq. Q. The dashed and dotted 
lines in each geometry are also conformally mapped to 
each other. 



Before presenting the desired conformal mapping, we 
discuss the implications of CI in these geometries. First, 
we present the implications of conformal invariance for 
the probability distribution of the domain wall in these 
two geometries. We can measure the probability pi(n) of 
the domain wall crossing the dashed line n times in the 
rectangle geometry, (a) , as well as the probability p[ (n) 
of the domain wall crossing the line n times in the slit ge- 
ometry, (d). In Fig.^ w e also show the mapping of a DW 
which crosses the dashed line three times. In addition, 
we can measure the probability p 2 (n) of the domain wall 
crossing the dotted line n times in the rectangle geome- 
try, as well as the probability p' 2 (n) of n crossings in the 
slit geometry. By statistical translational invariance, this 
probability is equal to pi(n). Thus, a naive application 
of CI would suggest that p^(n) — Pi(n). However, in the 
continuum limit it is not possible to distinguish between 
n and n + 2 crossings of the dashed line for n ^ 0; if 
the distance between successive crossings of the dashed 
line in Fig. ^ is of the order of the lattice spacing, this 
looks at large scale like a single crossing of the dashed 
line. Thus, the predictions of CI are: 
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(a) 



(d) 





FIG. 1: Mapping between different geometries, and of the 
DW, shown as the dashed-dotted line. 
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These identities in Eq. JQ) are the consequence of a con- 
formal automorphism of the slit geometry which inter- 
changes the dotted and dashed lines: z — > F(F~ 1 (z)+ir). 
We explicitly checked for a small rectangular system 
(N = 30*30) that pi — p2 ~ 0.603 and they are equal to 
p' 2 (quoted in tables 1 and 2) within the statistical error. 

Our main test for CI is to find whether p\ = p' 2 . The 
naive expectation is that the DW is less likely to cross the 
dashed line, due to the constriction reducing the number 
of possible configurations, while CI instead requires these 
probabilities to be equal. 

Constructing the Conformal Transformation — We de- 
fine the rectangle (a) to have width 2ir and height h, and 
thus aspect ratio h/2n. The horizontal direction as plot- 
ted is the real coordinate running from — n to n while 
the vertical is the imaginary coordinate, which runs from 
—ih/2 to ih/2. We define r(/i) = 2-Ki/h, and A(r) to 
be the modular lambda elliptic function. We define a 
function g(z) to map from (a) to (b) in Fig.^by 



g h (z) = sn(2izK(X(T(h)))/h\X(T(h))), 



(2) 



where K is the complete elliptic integral of the first kind 
and sn is the Jacobi elliptic function 8] . This function g 
maps the dashed and dotted lines as shown. The dotted 
lines extend off to infinity and all lines lie on the real axis. 
The upper and lower lines of the rectangle are mapped 
to the two solid lines, which have endpoints at ±1 and 

±i/vxRftjy. 

The mapping from (b) to (c) is simply z — > sz, for some 
parameter < s < 1; smaller values of s produces deeper 
cuts into the rectangle in (d) . In (c) , the endpoints of the 
lines are at ±s and ±s/a/A(t). We then determine the 
height h! of the rectangle in (d) such that 



l/V\^rW) = s/VWm- (3) 

The final mapping from (c) to (d) is z — > g^, (z). This 
maps the portion of the solid lines in (c) between ±1 and 



±l/y/X(rQi')) onto the upper and lower borders of the 
rectangle in (d) , while the portion of the solid lines in the 
third geometry between ±s and ±1 are mapped onto the 
slits. Thus, the full mapping from (a) to (d) is 



F(z)=g^(sg h (z)), 



(4) 



and the endpoints of the slits in (d) are located at 
F(±ih/2). 

We proceed by first finding the ground state of the sys- 
tem, using a mapping to a graph-theoretical problem, the 
minimum- weight perfect matching problem Domain 
walls were created by flipping the signs of the horizon- 
tal bonds in a column. Because of the periodicity in the 
horizontal direction, this induces a DW across the system 
and it is the crossings of the dashed central line and the 
dotted 'end' line which we study. A DW is best defined 
as a walk on the lattice dual to the original lattice, and 
the dotted and dashed lines are also lines of this dual 
lattice. 

In Tables 1 and 2 are displayed the values of p[ and p' 2 
and the probabilities of zero crossings, Pi(Q) and p' 2 (0) 
for s = 0.95 and s = 0.90 for various 'sizes'. Thus 30 * 32 
means that the system studied is rectangular with 30 
spins on each horizontal line (the direction in which the 
system is periodic) and 32 spins on each vertical line. 
The next number 8 is the number of rows cut by slits (4 
at the top of the system, 4 at the bottom of the system) 
indicating that there are (32 — 8) rows not cut by slits. 
(For small sizes it is not possible to find integers to get 
the sizes precisely correct for the given the aspect ratio). 

TABLE 1: Approach to CI for s = 0.95 
SIZE pi p' 2 pi(0) p' 2 (0) samples 



30 * 32; 8 
52 * 55; 14 
76 * 80; 20 
98 * 103; 26 



0.624 0.603 0.463 0.409 20000 

0.618 0.609 0.424 0.392 10000 

0.616 0.615 0.417 0.385 10000 

0.617 0.615 0.414 0.385 7000 



TABLE 2: Approach to CI for s = 0.90 



SIZE 


P'l 


P'2 


Pi(0) P2(0) 


samples 


30*33; 12 


0.639 


0.604 


0.494 0.407 


20000 


50 *56; 20 


0.623 


0.612 


0.448 0.399 


10000 


98 * 109; 38 


0.620 


0.611 


0.421 0.384 


6000 


124 * 138; 48 


0.612 


0.610 


0.411 0.381 


6000 



For each probability p in these tables y/(p(l— p)/N s ) is 
its standard deviation, where N s is the number of samples 
i.e. bond realizations averaged over. As the size increases 
so the continuum limit is approached, the closer p' x and 
p' 2 become, implying that the distribution of the DWs is 
conformal. In Fig. wc have plotted p' x versus 1 /L^ 1 
and p' 2 versus 1 / L 2 f 1 . (We have no proof that this is the 
way that p' x and p' 2 approach their asymptotic limit, but 
this dependence will be partly motivated below). Again 
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FIG. 2: Approach of p[ and p' 2 to asymptopia. Main figure is 
for s = 0.95 while inset for s — 0.90. 



one can see that in the continuum limit CI seems to hold. 
L\ is the number of rows not cut by the slit and L 2 is the 
total number of rows in the system. Thus for the 30*32 
system size, L\ = 24 and L 2 = 32. 

Unfortunately the probabilities for zero crossings p[(0) 
and p' 2 (0) appear to approach each other very slowly (if 
at all). In an attempt to understand this behaviour, we 
have studied the probabilities of n crossings, p[(n), of the 
central (dashed) line, and p 2 (n) of the right (dotted) line. 
p[(n) has a scaling dependence on the number of cross- 
ings n as l/-^]/ fa{n/Li ) and similarly p' 2 (n) is of 
the form l/L 2 f 1 f a {n/L 2 t The subscript a is added 
(so a — e or a — o) to allow us to distinguish even values 
of n from odd values, as the DW is topologically very dif- 
ferent depending on the parity of n. p' x [n ^ 0) can mean 
that macroscopically the DW crosses the central line once 
(say) , but then if one zooms in on that single crossing, it 
actually crosses many times. To see how many times it 
would cross, suppose the DW has fractal dimension df. 
Then, the intersection of the DW and a vertical line has 
dimension dt — 1. Thus, the expected number of cross- 
ings would be 1 which gives the above scaling forms. 
In Fig. |3 we display the data. The scaling functions for 
even and odd seem to be very similar, even though the 
DWs are topologically very different. Notice that the 
no-crossing probabilities are clearly not part of this scal- 
ing form (which they cannot be if they are non-zero). 
However, we suspect in the light of the above that the 
convergence of p\ (0) and p' 2 (0) to each other might be as 
slow as 1/L df ~ 1 : and as df is about 1.27 this could be 
a very slow convergence rate. We have also studied the 
case of s = 0.85, which corresponds to a still deeper cut 
by the slits. For this case, the convergence of even p[ to 
p' 2 has not been achieved in the largest sizes we have stud- 
ied (146 * 175), and we suspect that the slow convergence 
here is of similar origin. The conformal invariance found 
in the DW distribution encouraged us to find out if the 
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FIG. 3: Scaling of probability of crossings for s — 0.95. 

domain walls were also SLE processes Q. Suppose the 
domain wall is the curve j(t) which begins at a point on 
the boundary of the upper half plane H . The half-plane 
H minus the curve ^(t) can be mapped back onto H by 
an analytic function gt{z) which is made unique by de- 
manding that gt(z) ~ z + 2t/z + 0(l/z 2 ) at infinity. The 
growing tip of the curve is mapped onto the real point 
£(t). The DW is an SLE process if £(i) is a Brownian 
walk whose elements have an independent Gaussian dis- 
tribution and (£(£) 2 ) = K t. The diffusion coefficient n is 
of prime importance as it is related to the central charge 
of the conformal field theory |7j. In practice we approx- 
imate gt(z) by composing a sequence of discrete, confor- 
mal slit maps of the form z — > y/(z — £i) 2 + 4Ati + £j , 
where the parameters Ati , Ci are chosen so that the i-th 
such map removes the i-th step from the domain wall 
following the procedures of Ref. 0, to produce a se- 
ries of times U = ti—i + Ati and values £(ti) = £j which 
approximate £(i). We denote the coordinates of the do- 
main wall points by zf, i = 1...N. The first such slit map 
transforms the coordinate of the first step, z®, into the 
origin, and transforms Zj, j > 1, into a new point z 1 . 



In general, the j-th map transforms 
and gives 



into the origin 



(5) 



-112 



(J > *)• 



As usual the complex number z = (x,y) and £(to) = 0. 
The sign of the square root is chosen so that it has the 



sign of [x % - 



The geometry which we studied 



was a square L x L with periodic boundary conditions in 
the horizontal direction and open boundaries in the ver- 
tical direction which is the direction in which the domain 
runs. In Fig. 0]we show the average over realizations of 
the disorder of (£(i) 2 ), plotted against t, for three values 
of L. When L = 180, we took 3000 disorder realizations, 
L = 220, 4000 and for L = 300 we took 5000. (£(t) 2 ) 
is linear in time for a range of times which increase with 
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the system size L and from the slope of this linear region 
we estimate that k ~ 2.1. Our boundary conditions, to- 
gether with the finite size of the system, means that it 
does not properly satisfy the requirements for produc- 
ing either chordal or dipolar SLE 7], which may partly 
explain the modest size of the linear regime. In a recent 
related study it was found that u sing the dipolar SLE did 
indeed extend the linear regime [llj. In the inset to Fig. 
Hlwe show that the probability distributions of £(/)/\/fti 
at four different times within the linear regime are stan- 
dard Gaussian as would be required for the domain walls 
to be SLE processes. 

If spin-glass domain walls are SLE processes, there may 
be a relationship between the fractal dimension of the 
DW df and the exponent 9. The fractal dimension is 
related to k via df = 1 + k/8 [7j. (Our numerical value 
for k ss 2.1 and the estimates of df ~ 1.27 ± 0.01 in Ref. 
are consistent with this relationship) . The correlation 
length exponent exponent v is related to one of the Kac 
elements of conformal field theory: for example in Potts 
models with components Q, 1 < Q < 4, one has d— \jv = 
2 + 9 = 2/12,1, but in other models d — 1/v is 2/11,3 or 
2/11,2- d is the dimensionality of the system i.e. 2. Now 
if SLE applies, each of these elements is related to n: 
2/i2,i = (6— k)/k 0|, 2/11,3 — k — 2 etc. For each of these 
possibilities one can derive a relationship between df and 
9 and the only one which comes close to the numerical 
values d f = 1.27 ± 0.01 and 9 = -0.285 ± 0.05 is from 
2/12,1 = — - = d — - = d + 9. Then on eliminating k in 
favour of df gives 



0.006 



4(3 + 6») 



(6) 



On using one of the alternative possible relations, say 
2 + 9 = 2/11,3, d f = (12 + 60/8. Then the predicted value 
of df = 1.46 which is not consistent with its numerical 
value. Eq. © seems to be the only possible relation- 
ship between df and 9 which is compatible with their 
numerically well-established values. (Note that Eq. 
would not apply to the ±J spin-glass model as for it v 
might have the same value as for the spin glass model 
whose bonds have a Gaussian distribution |l2j but it has 
9 = 13]). The apparent success of Eq. EJin provid- 
ing a relationship between df and 9 might provide a clue 
in finding the kind of conformal field theory appropriate 
for two-dimensional spin glasses. Our numerical evidence 
for conformal invariance and SLE strongly suggests that 
such a field theory should indeed exist. 
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FIG. 4: (£(t) 2 ) versus Loewner time t for three different sys- 
tem sizes L. The straight line has slope 2.1, which is our es- 
timate for k. Inset: probability distribution at four different 
times, 0.0003, 0.0004, 0.0005 and 0.0006 when L = 300. The 
curve is the Gaussian exp(— x 2 )/\/2tt, where x = £(t)/\/fiJ. 
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